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الباب الخامس الملاءمة والانكفاء بواسطة البرمجة 
5 -1 الملاءمة و الانكفاء OO A i a a a a a a A a‏ 


كلمة مدير الموقع 


بسم الله الرحمن الرحيم 


الحمدلله والصلاة والسلام على أشرف الخلق والمرسلين سيدنا محمد صلى الله 
عليه وسلم .. 
إن المتابع الجيد للساحة العربية على الشبكة العنكبوتية ليجد أنها تفتقر الى 
الكثير من المواضيع التي تهم المتخصصين كالمهندسين والفنيين والتقنيين 
عموماًء وقد نجد في بعض المواقع بعض الاجتهادات من توفير المراجع ولكنها 
في القالب كون مصار اك وبلغات قير اللغة العرية. 


من هذا الباب .. عكف فريق موقع و منتديات التقنية ( أكبر تجمع للمهندسين 

العرب) أ6١.6٣‏ ).ر على تبني مشروع بناء مكتبة عربية متخصصة في 

المواضيع الهندسية والتقنية › بهدف نشر هذه الكتب وتوزيعها الكترونياً عبر 

الموقع لتساهم ولو بجزء بسيط في اثراء الساحة العربية بالمواضيع العربية 
الهندسية. 


وهذا الكتاب الالكتروني الذي بين يديك هو نتاج هذا المشروع الجميل › الذي نسأل 

الله أن يمدنا بالعون والتوفيق على اكمالهء والشكر لجميع الاخوة ممن ساهم في 

اخراج هذا الكتاب أو ساعدنا بالرأي والمشورة والاقتراح وصلى الله وسلم على 
سیدنا محمد وعلی آله وصحبه وسلم . 


بقلم المهندس/ فهد الرفاعي _ 
أستاذ في كلية جازان السعودية 

والمشرف العام لموقع التقنية 
refaeefa@hotmail.com‏ 


الباب الأول 


مقدمة في التحليل العددي 
1 -1 استهلال introduction‏ 


التسمية الدقيقة لهذا التخصص هو التحليل العددي( numeric]‏ 
ئ«اهمه) يستخدم التحليل العددي في حل المعادلات الرياضية التي يصعب حلها 
أو يستلزم وقت طويلا في الحل» فإذا كانت لدينا معادلة رياضية من الدرجة 
الرابعة أو الخامسة فان من الممكن حلها بالطرق التقليديةء و حتى هذه الطرق 
تاخد وقتاء إما إذا كانت المعادلة من الدرجة السادسة و أكثر يكون من الصعب 
اتضحت أهميته البالغة في حل هكذا معادلات» و ذلك لتوفيره الوقت و الجهد. و 
بالأخص في المعادلات التي تحتاج إلى تكرار كبير من اجل الوصول إلى النتيجة 
أو الحل. 


عندما يحصل فني أو مهندس على وظيفة معينةء مثلا في الصيانة أو 
الإنتاج أو الجودة أو تدبير المواد الأوليةء يتوقع آنه سيفتح جهاز الحاسوب ليجد 
أمام4 İحد yi Oracle gi PeopleSoft yi Baan İi SAP al‏ 
۴65ل و كل ما عليه هو الضغط على بضعة أزرار هنا و هناك ليعرف 
حالة المخزون و الطلبيات و موعد تسليمها و نتائج اختبارات الجودة و كمية 
المواد الأولية التي يجب شراؤها و عدد ال ءاه| في خطوط الإنتاج و أسماء 
العمال الذين عملوا في تجميع منتج معين و عدد ساعات غياب عامل ماو 
الاقتراحات التي قدمها العمال و عدد توقفات آلة ما ...لكن هذا غير موجود على 
أرض الواقع» واقعنا العربي... 


و حتى الشركات الغنية كشركات البترول و الصلب» التي استثمرت في هذه 
الحلول و ركبتهاء لا تستفيد منها كاملا و أغلب الموظفين يعتبرونها عبنا إضافيا 
على واجباتهم اليومية و قد يستعملونها فقط عند قدوم زوار خارجيين. أما في 
الشركات المتوسطة و الصغيرةء فمن الصعب على الفني أو المهندس إقناع 
المدراء بالفائدة الملموسة لمثل هذه الأدوات خصوصا أن ثمن تملكها السنوي 
Cost o Ownership‏ اtaدآ‏ مرتفع مقارنة بميزانية الشركة السنويةء لذلك 


فالحل الوحيد الذي يجده المهندس لتدبير أموره اليومية هو ابتكار حلول برمجية 


و لن يتأتى له ذلك إلا بمعرفة مسبقة بطرق التحليل العددي و استعمال 
الحاسوب في حل معادلات رياضية. و هذا يدرس في الجامعات و مراكز التكوين 
(و منتدانا ١‏ لحبيب) و ليس في مكان العمل. لان وقت العمل ضيق و محسوب .كما 
أنه لا مجال لاستعمال برامچ مقرصنة تحت طائلة تعريض الشركة لغرامات مالية 
في حال قيام مايكروسوفت و أخواتها بالتحقق من برمجيات الشركة (إثر مكالمة 
هاتفية من عامل لم يستفد من العلاوة السنوية ). 


E E E E E E EO 
تستعملها ۀ منزل لا تس نثبیتھ كلها في حاسوب العمل دون‎ 


في النهايةء أقول إني لا أقصد أن يكون المهندس مبرمجا محترفا زيادة 
على التخصص الذي درسه»ء لكن وضع برنامج صغير ل نبع مؤشر ما او ورقة 
امء×ه ببعض المعادلات الرياضية أسرع و أنجع من تعلم برامج من العيار الكبير 
في حالات كذيرة. 


في حقيقة الأمر يوجد الكثير من البرامج يمكن استخدامها في هذا المجال 

مٿل برnli MATHCAD , MATLAB‏ و ت ڊرgnli EXCEL‏ 
المرlفق‏ نحajة OFFECE‏ ن أجر اءه في حل الكثير من المعدلات و إجراء 

العديد من الصغير الرياضية المفيدة › لكن في هذه السلسلة سنقوم باستخدام إحدى 
لغات البرمجة في التعامل مع المعادلات الرياضية > سنتعامل مع لغفة 
۴ORTRAN‏ لما تمتاز به هذه اللغة من استقرار عال و دقة ق المعادلات و 

تحديد دقيق لنوع المتغيرات مع سعة كبيرة في نوع ا رت من الممكن أنها 


ليست على درجة كبيرة من الانتشار و الاستخدام › > غير أن اأ تخصصين في مجال 
لتحليل العددي و الرياضيات بصفة عامة يعترفون لها بالفضل» سنستخدم الإصدار 
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و هو إصدار يعمل تحت بيئة الويندوز(۷]N(0WS‏ ) و لست في معرض 
الدعاية لهذه اللغفة و لاأحب المقارنة بين اللغات المختلفة غير أن 


۴ORTRAN LANGUAGE‏ هي لغة المهندسین بکل کفاءة. 


غير أن هذا الكتاب لن يهمل بعض اللغات ذات الانتشار a‏ 
و قد اعتمدت بعض الطرق المدرجة في هذا الكتاب بلغة ٥‏ من اجل خلق تنوع يفيد 


قبل التعامل مع المعدلات الريات 
التواب و المتغيرات. 


FOR T8۸ Nنارتروف أساسیات لغة‎ 


CHARACTERS jgaرll‎ 1 


تستخدم مجموعة من الرموز الأساس 


في لغ ۴ORTRAN‏ و تتکون 


4 


\CTERS مlقرألا‎ 1.1 


RS CHAR 


NUMB 
9 8 7 6 5 4 3 2 1 0 وهي‎ 


" 


أنواع البيانات في غ F۴ORTRAN‏ 


تمثل البيانات لفي لغة الفورتران ۴0۸۲۸4۸۸ بإحدى الأساليب الآتية : 


RL INFORNATION- e المعلومات‎ 
› العمليات الحسابية‎ a | م‎ ek تست ثيل‎ 
(TRUE OR FULSE) ° 


1-2 طيقة المقاطع 
BISECTION METHOD‏ 
لنفرض أن لدينا المعالة الرياضية التالية : 
X^2 +4 =0‏ 


مع ملاحظة أن الرمز(^) يعني أس أي مربع القيمة ( ×) و 5^× يعني الأس 


سيكون الحل المثالي لهذه المعادلة هو كالأتي: 
X۸24‏ 
X=+14‏ 
2 
2-=) وهو حل دقيق تماما بمجرد إتباع الخطوات السابقة, نحصل على الحل. 
و إذا كانت لدينا المعادلة الرياضية التالية : 
X^3+4X^2-10=0‏ 


فحلهايرهق قليلا و لا يعطي قيم دقيقة › و لهذا نلجأ إلى طريقة المقاطع 
B [ECT] 0N METHOD‏ التي تتلخص فيما يأتي : 


لنفرض أن لدينا دالة (»)۴ معرفة و مستمرة خلال الفترة إ 8 ,۸ ]) مع 
الدالتین (۴)۸ , (۴)8 فهنالك قيمة 0× تنتمي إلى الفترة 


}۸,8{ كالدالة ۴))0(=0.0 » و التي تکون حلا للمعادلة ())۴ 
بما أن القيمتين 8 ۸ تنتمي إلى نفس الفترة فهذا يعني ان الدالتين 


(۴)۵ , (۴)8 تنتمي إلى نفس النطاق و يكون الرسم البياني للدالة كالآتي. 


FPIONT SOL UTION 


(شکل2 -1) 
و حل المعادلة يكون هو النقطة 0× ( P]ON† S0107110‏ ) و بما أن 
أساس الحل تقريبي فتكون الطريقة المثالية لخطوات الحل هي اخذ مجموعة من 
النقاط التقريبية المقربة للقيمة الصفرية كالقيمة 0.00001 
و اختبارها فهي لا تساوي صفر لكنها دقيقة جدا و مقربة جدا من الصفر حيث 
يصعب الحصول على الحل الصفري في التحليل العددي مباشرة و هذه القيمة 
سنفرضها ليحملها المتغير ۴ و يقترب كل مرة بقيمة دقيقة من الحل الصفري و 
تعتمد طريقة المقاطع M٤1٤1 0D(‏ 10۸Nآ)٣B]SE)‏ على نفس الفكرة في 
التعامل مع المسائل الرياضية و خوارزمية البرنامج هي كالآتي : 
START‏ 
RAED (A,B)‏ 
DO P=(A +B )/2‏ 
Y1= A^3 + 4*A^2 - 10‏ 
Y2= B^3 + 4*B^2 -10‏ 


YO0=P^3 + 4*P^2 - 10 


IF (¥0*Y1 > 0.0) GOTO STEP 10 
A=A 

GOTO STEP 20 

DO A=-P 

B=B 

GOTO STEP 20 


WRITE Y¥0,P 


الاختبار الأول يسال هل قيمة 0.00001 >۲0 بم 
ال اډ 4 ذعد ۳ ل al‏ الاخ 8 ا 


SOLUTION C THIS PROGRAM TO CALCOLATE 
METHOD OF EQUATION BY USED BISECTION 


P=(A +B )/2 

Y1= A^3 + 4*A^2 - 10 

Y2= B^3 + 4*B^2 -10 

Y0=P^3 + 4*P^2 -10 
IF ( ASB( Y0) .LT. 0.00001 ) GOTO 100 
IF (YO*Y1.GT. 0.0) GOTO 10 
A=A 
B=B 
B=B 


RIT] E) Y0 ,P 


Nuten neth0d نjننروین ر طريفَة‎ 


في هذا الدرس نستعرض سويا احد الطرق المتبعة في التحليل العددي بعد 
أن تكلمنا عن مفهوم التحليل العددي و أهميته و الطريقة الثانية المتبعة في 
التحليلي العددي هي طريقة نيوتن (04ط٤٤ہ ٤٥١‏ سم ) و لا أود أن أتطرق 
إن الاستتتاج الرياضن لهت الطريقة يالقاس الذي ارغب فيه في الترکیز على 
البرامج الرياضية لها تعتبر هذه الطريقة من الطرق السهلة في إيجاد القيم 
التقريبية للمعدلات الرياضية و الاختلاف الأساسي بيناه و بين الطريقة السابقة ( 


طريقة المقاطع) (ك0طاء. «٠ااءعءط)‏ هو سهولة استخدام طريقة نيوتن و 
سرعتها في إيجاد الحل التقريبي للمعادلة. 


المفهوم الرياضي لنيوتن 


اذا كانت لدينا دالة حقيقية (»)۴ و كان 0× هو الجدر للمعادلة المطلوب الحصول 
عليه ( الحل التقريبي للمعادلة) فيمكن إيجاده عن طريق الآتي: 

X2=x1*f(x)/f '(x) 

حیث : 

2× هي القيمة التي نبحث عنها أي نفس قيمة الجدر ( قيمة 0× ) 

1× هي القيمة المدخلة عند القراء 

(×)۴ هي قيمة المعادلة بعد تعويض بي قيمة 1× في المعادلة 

(×) ' ۴ هي قيمة المشتقة الأولى للمعادلة بعد. 

و تكتب خوارزمية البرنامج كالأتي: 


Read x1 1 
f(x1)=x1^3+4*x1^2-10 .2 


f ' (x1) = 3*x1^2 +8*x1 .3 

if f(x1) > 0.00001 write x2,f(x1) .4 

x2=x1-( f(x1)/ f ' (x1) ,x1=x2 goto .5 
step (2) 

write x2,p .6 

stop .7 


شرح الخوارزمية 

يبدأ البرنامج بقراءة قيمة الدالة عند النقطة (1× ) ثم تتم عملية تعويض من اجل 
إيجاد القيمة الفعلية للدالة عند نفس النقطةء تعوض قيمة (1× ) في مشتقة الدالة و 
كأننا أوجدنا الميل في هذه الحالةء تتم اختبار قيمة الدالة (1») ۴ فإذا كانت أصغر 
من 0.00001 فتتم طباعة القيمة مباشرة و إلاننتج قيمة جديدة هي 2× من 
طرح قيمة 1× من مقسوم الدالتين (×) ۴, (× ) ' ۴و هي قيمة أو نقطة تقاطع 
المماس مع محور السينات و عندما نجعل قيمة 1=×2× فإننا نقترب من الحل أكتر 
أي من النقطة الصفرية التي تحقق الحل. 


و هذا هو البرنامج بلغ -:FORTRAN‏ 


Read(*,*) x1 

10 f(x1)=x1^3+4*x1^2-0 
f ' (x1) = 3*x1^2 +8*x1 
if f(x1) > 0.00001) goto 20 
x2=x1-(f(x)/f ' (x)) 
x1=x2 

goto 10 

20 Write(*,*) x2 ,p 
stop 
end 


و يمكن استخدام جملة 00 بدلا من جملة 6010 لأداء نفس المهمة على 
النحو الآتي: 


Read(*,*) x1] °“ 
DO 500,I=1,N ° 
f(x1)=x1^3+4*x1^2-10 


f ' (x1) = 3*x1^2 +8*x] 

if f(x1) > 0.00001) GOTO 20 ° 
x2=x1-(f(K)/f ' (x)) ® 

x]=x2 ® 

GOTO 500 ° 

500 CONTINUE ° 
20 Write(*,*) x2,p ° 
stop ® 

end ° 


و الرسم التالي يوضح فكرة نيوتن في حل المعدلات : 


و نلاحظ انه بدأ من نقطة بعيدة و في كل مرة يقترب بمقدار معين من القيمة 
المطلوبة 0> و كأنه يأخذ الشكل التالي: 


(شکل 3-2) 


2 -3 طريقة نيوتن رافسن ونيوتن رافسن المعدلة 


Newton-Raphson Method 
And 
Modified Newton-Raphson Method 


طريقتا نيوتن رافسن ونيوتن رافسن المعدلة: 
هما طريقتان لحل المعادلات الآنية غير الخطية , ولو أخذنا معادلتين في 
متغيرين كمثال ونريد ايجاد قيمتي رر× اللتان تحققان المعادلة (1) 
F1(x,y)=0 F2(xy)=0............ (1)‏ 

أولا: طريقة نيوتن رافسن 
من متسلسلة تايلور في متغيرين كمثال 

FIGYuD=FI(xpyD+ (u xD+ (yp)... 

F2 D=F2(xpyD+ E (ua x)+ (pi)... 


نفترض أن 1+1ر,1+× فريبة جدا من الحل لذلك 
F1(xi+1,yi+1)=F2(xi+1,yi+1)=0‏ 


و صم كذلك 
X1 = Xi =h‏ 
Ji1” yi= k‏ 
لحصلنا على الصورة 
=F1(x;y;+ E h+ Fk‏ 0 
0x dy‏ 
=F2(x;y;+ 2 h+ 2‏ 0 
0x dy‏ 
ومنها 


۴1 p+ °F ( = - F1(yy) 
0x dy 


F2 ْک‎ k = - F2(xy;) 
ر‎ 


0x 
وهذه المعادلة غير متجانسة ونحصل على حل عندما محدد ١141طc0 4ل لايساوي‎ 
الصفر‎ 
OF dF 
e e) 
j= F24 F24 0 
dy dy 
ونحصل على المعاملات كالتالي:‎ 
- FI(i, yî) u 
) 
H= J 
- F2(i, yi) e 
OF. =F) 
K= 3 
F4. = F2(i yi) 
نستخدم المعادلتين التاليتين للحصول على قيمتين جديدتين لكل من ,× غير‎ 
الافثر اشيشن‎ 
Xıı=xi+h 
Yi,ı=yi+k 


نعيد الكرة من جديد باستخدام القيمتين الجديدتين . 

ثانيا: طريقة نيوتن رافسن المعدلة 

من الواضح صعوبة الطريقة السابقة في حالة تزايد عدد المعادلات (م) أي 
سيكون لدينا (”«) من المشتقات الجزئية لكل معادلة «” من المشتقتات 
والمتغيرات , لذلك تعدل الطريقة كالآتي: نحصل على قيمة × الجديدة من طرح 
الحالية من إحدى الدالتين على تفاضلها الجزئي بالنسبة ل × أما قيمة ر الجديدة 
فمن طرح الحالية من الدالة الأخرى على تفاضلها الجزئي بالنسبة ل لإ . 
واختيار أي الدالتين نستخدمه مع × أو ر ليس عشوائيا إنما يكون الاختيار 


مهما حيث يؤدي الاختيار الصحيح إلى تقارب نحو الحل أما الخاطئ فسوف 
يبعدك عن الحل. وهنا کمتال اخترنا ۳1 مع ×و ۴2 مع ل 
F1(xi, yi)‏ 
)1( و 
dx‏ 


ED 
Yi+1=Y¥i OF 1C, yi) eeeececceececcccccccecccecccccccscsese )2( 


dx 


Xi+1=Xj” 


بالنسبة للدالتين : (×)«iو+×=ر‏ و 0=ر-(3si")x‏ 
واللتان يراد نقاط الحل لهما أي نقاط تقاطعهما ويبين الرسم في الصفحة التالية 
نقاط التقاطع. وسنقوم بايجاد الحل بالطريقتين: 


بطريقة نيوتن رافسن 
F1(x,y)=x+sin(x)-y‏ 
F2(x,y)=3sin(x)-¥y‏ 
نفاضل كلا الدالتين بالنسبة ل× وبالنسبة لر 
ات -1+C05(x)‏ ۴1 
dx dy‏ 
ت (×)3¢05= °۴2 
0x dy‏ 
1+cos(@) -1‏ 
J= = -(1+cos(x))+3cO0s(x)‏ 
3cOS@) -1‏ 
SET EONKE TS. sois (i‏ 
(x + sSin(x) — y) -1‏ - 
H= /J = [(x+sin(x)-y)-(3sin(x)-y)]/J‏ 
1- ( ر = (3s1n()X)‏ - 
E OL (i)‏ 
cOS(x) - (x + sSin(x) ¬ y(‏ +1 
K= /3J =‏ 


308) ¬ (3s1n(X) — y) 


5 EERE 4F 


=) E -2yc0ئ‎ ey 
= ( y-3sin(x)+cos(x)( 3x-2y ))/J .. Sonne U) 


1X, y 
تمثل النقطة الجديدة‎ :Xnew, رnew‎ 

عامل جاګويي وی اوي المعادلة (1) 
4 معامل تفاضل الدالة بالنسبة 


k=(y-3*sin(x)+cos(x)(3*x-2*y))/jac g h=(x-2*sin(x))/jac Jaج‎ 
ynew=y+k gy xnew=x+h Jaجl‎ 

9. أضف 1 إلى 1 

)1. اطبع قيم لإ ر× , 1 الجديدة 


1. هل القيمة المطلقة للفرق بين × الحالية والجديدة أقل من ٤‏ إذا كان 
لا..... استبدل × الحالية بالجديدة انتقل إلى 5 
2. هل القيمة المطلقة للفرق بين ر الحالية والجديدة أقل من ٤‏ إذا كان 


لاء.... استبدل ر الحالية بالجديدة انتقل إلى 5 
3. النهاية 
المخطط الانسيابي 


Flow Chart 


X=xnew 


= , 
Tar—2*rcac(v)-1 
yes 

NO : 


H=(x-2*sin(x))/jac 


K=(y-3*sin(s)- 


Divide by zero 
Xnew=x+h Tacrnhin—0 


Vnow—vkk 


no e 
yes 
no n 


T 
yes 


# include<conio.h> 
# include<stdio.h> 
# include<math. h> 
void main () 

1 LE *stream; 


ynew, jac, h,k; 
ntf ( initial value of x, y :"); 
1f("%Sf Sf", &xnew, &ynew) ; 


NF 
stream = fopen ("newraf.FIL", "w+"); 
fprintf (stream, "The initial value of x,y is (%2.1f 
«2.1 E 7) Tas 
RE ê sam, "\n i\lt x\t\t y\t"); 
KT OE i. "7; 


fprintf (stream, "\n %d\t %f\t Sf\t", i, xnew, ynew) ; 


if( jae ==0) { 
fprintf (stream, " divide by zero , Jacobin=0 "); 
aE و‎ 


fclose n 2 


أن هناك ثلاثة جذور أحدهم (0,0) والآخران أحدهما موجب 


سالب لڪل من ز× ومن هنا 


Enter initial value of x, y :2 1 
The initial value of x,y is (2.0,1.0) 


2. 000000 


000000 
e 


: 843207 
. 843241 


` WN Ho 
bhhHhHH N 


494 895 و1 
Enter initial value of x, y :-2 -1‏ 


The initial value of x,y is (-2.0,-1.0) 


-2. 000000 -1. 0 


E -2 47 
895494 -2. 1 
-1.895494 -2.441 


ھە دا دا س ک 


Enter initial value of x, y :0.5 0.5 


The intial value of x,y is (0.5,0.5) 


-0.000000 -0.000000 
0. 000000 0. 000000 


` Ww N HH O | 


أما الطريقة المعدلة 
لو أخذنا الدالتين بهذا الترتيب 


F1 (x,y)=3sin(x)-y 
F2 (x,y)=x+sin(x)-y 


وال ناراي قال باشبة ل وة يبةن 


F2 1 


F1 . 
dy 


=3c05)%( 


X;, = xi - (3sin(xD-yi)/ (3cos(xD) 
Y;,ı = yİ - (xi+sin(xi)-yi)/(-1) 


.10 


0 VO 5Û ضۅ‎ i b0 ا‎ 


٠ أبد‎ 

ادخل قيمة ر , × الأولية. 

اجعل 1=0 وهي تمثل عدد التكرارات 

اطبع قيمة ر ,× ,1 الأولية. 

ynew=y-(x+sin(x)-¥y)/(-1) gy xnew=x-(3*sin(x)-¥)/(3*c0sS(x)) Jعجl‎ 
1 أأضف 1 إلى‎ 

اطبع قیم ۷W٥٣ر‏ ,Wعہ×‏ ,1 التي تمثل قیم ,× الجديدة 

هل القيمة المطلقة للفرق بين × الحالية والجديدة أقل من ٤‏ إذا كان 
لا..... استبدل × الحالية بالجديدة (Wع«»)‏ ثم انتقل إلى 5 

هل القيمة المطلقة للفرق بين « الحالية والجديدة أقل من ٤‏ إذا كان 
لاء... استبدل ر الحالية بالجديدة ٣٥ ۷W(‏ ر) ثم انتقل إلى 5 


المخطط الانسيابي 
flow chart‏ 
| 1 1 


X=xnew , 


Xnew=x-(3*sin(x)-y)/(3*CcOS(X)) 


Ynew=y-(x+sin(x)-¥y)/(-1) 


و 
e‏ 


T—-Tı1 


n 


yes 
no 
yes 


End 
)5- 2 (شکل‎ 


# include<stdio. h> 
# include<math. h> 
void main () 
{ FILE *stream; 
float x, y, xnew, ynew; 
printf ("Enter initial value of x, y"); 
scanf ("%f Sf", &xnew, &ynew) ; 
1=0; 
fprintf (stream, "The initial value of x,y is 
($2.1f ,%2.1f£) \n", xnew, ynew) ; 
fprintf (stream, "\n i\t x\t\t y\t"); 
fprintf(stream,''ln_________ ); 


Sf", i, xnew, ynew) ; 
{ x=xnew; 
y=ynew; 
i++; 
Sf", i, xnew, ynew) ; 


fclose (stream) ; 


} 


Enter initial value of x, y 


The initial value of x,y is (2.0,1.0) 


e 
2.735107 


. 881608 
۰ 80 E 


2 
.86343 
. 81 


. 8 
. 2 


0وو 
E‏ 


2 


738 
870 


7و ` 
880843 . 
0 . 


85076 
. 4 
. 0 
.836 
.847000 
5 


4م 
8457. 


.841903 
. 7 
.844199 
. 844 
.842563 
.842748 


.8428 
. 2 
.84346 
.843419 
. 847 
.845 
. 845 
. 843331 
.843153 
. 8437 
. 844 
. 847 
e 
` 8423 
 .N 4 
(1- 2 (جدول‎ 


E N O INNS NNO 


. 906570 
. 903323 
.887941 
. 884 
. 90106 
e 


03 
kk 
. 894064 


. 1 


. 896511 
.89499 
. 0 
. 896206 
. 896009 
. 890 
. 895128 
5 
. 5 
E 
.895309 
. 895676 
. 895626 
e 


54 


1 
¬: 7 
. 1 
. 8 
. 895461 


05 0.5 


ّ 


Enter initial value of x, 


The initial value of x,y is (0.5,0.5) 


0.01120 


y 


. 50000 


|. 3 
. 65 
.--.- 8 


. 5 


7 


ا 


07 
.014556 


. 1 
. 01 


. 0251 


. 0 0 


.002874 
. 00 
.001916 
.005013 


. 007 
.0032 


. 006466 


CSL e 


× 


0. 500000 
0.32886 
0.217908 
0. e 


. 042838 
. 0124 


.007278 
.01906 
. 004850 
Ol 


., 5 
. 00560 
0.001437 


O O 


. 003760 
. 000958 
. 0 057 
0.000639 


Oo © O 


0.001671 
0.000426 


3 


ONU WN HH © 


30 0.001114 
31 0.000284 
32 3 


E 


5 
.000378 
. 0000 
. 000252 
. 000660 
. 000168 
. 000440 
. 000112 
. 000293 
.00075 
«0198 


000330 .0 36 
0.000084 37 
0.000220 38 
0.000056 39 
0.000147 40 
0.000037 41 
0.000098 42 
0.000025 43 
0.000065 44 
0.000017 45 
0.000043 46 
0.000011 47 
ا 48 


yS 


. 000033 
OT 


. 000058 
. 000015 
. 000039 
. 000010 


51 0.000005 
52 0.000013 


e 


)2-2 (جدل‎ 
Enter initial value of x, y :-2 -1 


The initial value of x,y is (-2.0,-1.0) 


1 × y 


0 -2.000000 -1.000000 
1 -3.384041 - 27 
2 -2.137745 - 1 
3 -1.757074 -2 9 
4 -1.697342 > 4 
5 -2.321275: - 26 


57 -1.895541 - 209 
58 -1.895528 2 3 


= xi — (3sin(x;)-yi )/ (3*cos(xi)) 
= yi - (xı ıtSin(&i,1)-¥yİ)/(-1) 


و 
ا 


xnew=x-— (3# : 
ynew=y- (xnewtsin (xnew) -y) /(-1); 


٠)) (ر-‎ 7 (3 *COS (x) 7ر‎ 7 


فكانت النتائج لنفس الإدخالات كالتالي: 
Enter initial value of x, y :2 1‏ 


The initial value of x,y is (2.0,1.0) 


000000 .1 000000 . 
3.143961 384041 . 
942 .2 0577 . 
2.81201 س 7 
9 1 ۰ 2.8 6 


"0 
. 0 .8 

2ق 2 0 . 
2.841 9657 . 
2.831 88909 . 
867 .2 955 . 
4 .2 83 . 
9 .2 5 . 
5 .2 896 . 


“I GON QQ ÛU WwW N HF © 


2.844028 1.896649 14 
5 .2 1.894678 15 
849 .2 1.896078 16 
2.842959 1.8950 17 
2 .2 11.80 18 
2.843099 1.85 19 
2.843 1.854 .20 
2.843169 1.895388 21 
2.843 0 22 
2.845 1.11 :23 
847 .2 112 24 
2.843 1.895467 25 
2.845 1.84 26 
842 .2 1.895481 27 
2.843248 4ق 1 28 
2.8427 1.895487 29 
2.845 1.895499 30 
(جدول 3-2) 


Enter initial value of x, y :0.5 0.5 


The initial value of x,y is (0.5,0.5() 


3 × y 

0 0. 500000 0. 500000 
1 0.143613 0.286733 
2 0.095576 0.191007 
3 0.063669 0. 1275 
4 0.042432 0.084850 
5 0.028283 0.056563 
6 0.018854 0.037708 
7 0.012569 0.025138 
8 0.008379 0.016759 
9 0.005586 0.011172 
10 0.0074 0.007448 
11 0.002483 0.004966 
12 0.0055 0.003310 
13 0.001103 0.002207 


001471 . 
000981 . 
000654 . 
000436 . 
000291 . 
000194 . 
000129 . 
000086 . 
000057 . 
000038 . 
(جدول2 -4) 
و بتغيير اختيار الدوال(إعادة الترتيب) كالاتي: 
F1(x,y)= x+sin(x)-¥y F2(x,y)= 3sin(x)-y‏ 
ونأخذ للأولى تفاضل بالنسبة ل × والثانية بالنسبة ل ل 
7 
dy‏ 


. 000736 
. 000490 
. 0007 
. 000218 
. 000145 
. 000097 
. 000065 
. 000043 
. 00029 
. 000019 


14 
15 
16 
17 
18 
19 
20 
21 
22 
23 


OOOOCOCOOCO0CO0O0OO 
O OOOOOO0OO 


arı 


=1+c0sS(x) 
dx 
ونطبق القانون‎ 
Xj = xİ — (xi+sin(xi)- yi)/ (1+cos(xi)) 
Y;,1 = yi - (3sin(xi)-yi)/(1) 
فإن الناتج سيكون مبتعدا عن الحل وهذا هو ما حصل فعلاً‎ 


0 3. 000000 2. 000000 
1 -111.0243 0.4230 


71 -1431034113556480. 000000 -1. 9 
2 49127493789024256. 00000 -0. 1 
73 +NAN +NAN 


الاستنتاج 

* في الطريقة المعدلة وبعد استخدام قيمة < الجديدة لإيجاد ر الجديدة تقلص 
عدد التكرار إلى النصف تقريبا عن الطريقة المعدلة التي استخدمت فيها قيمة ب 
القديمة لإيجاد ر الجديدة أي وجدنا الحل في نصف الوقت. 

*في طريقة نيوتن رافسون تم الوصول إلى الحل بأقل تكرار من المعدلة 
وبتقارب دائم نحو الحل بينما الطريقة المعدلة نجدها مرة تتقارب وأخرى تبتعد 
مما يزيد عدد التكرار. وهذا يتضح من الأشكال التوضيحية لاحقا. 

* اختيار أي معادلة نعتبرها ((»)۴1 مهم. حيث نلاحظ التقارب والتباعد الناتج 
عن ذلك 

* نلاحظ من الجدول, اختلاف نتائج الطرق سالفة الذكر من حيث عدد التكرارات 
حسب قيمة رر× المدخلة مع ملاحظة أي من الحلول الثلاثة نصل إليه في كل 
مرة. ومن المعلوم أن هذه التكرارات تعتمد على 5 (قيمة الخطأ المسموح به) . 


(جدول 2 -5) 


نيوتن رافسن المعدلة مع نيوتن رافسن | , , . 
1 نيوتن راقسن ٍ 
استغلال × الجديدة لإيجاد ر المعدله النقطة المدخلة 
التكرار < الحل التكرار < الحل | التكرار < الحل 
24 < )0,0( و 5-,1.8-) | 5-,3- 
4 < (0,0) 54 < (0,0) 5 > (0,0) 0.6,0.6 
-,1.8-( 
 )1.8,2.8(<5 2.8(<5 (-1.8,-2.8(<33‏ 1.5,0.75 
0,0(<27( تباعد  )1.8,2.8(<5‏ 1.5,2.5 
P04 (1.8,2.8()>4 ٠ (1.8,2.8(<7 (1.8,2.8(<29‏ 
 (1.8,2.8(<61 (0,0(<24‏ 1.8,2.8(>5( 3,5 
0,0(<46( تباعد ‏ /1.8,2.8(<114) 7,4 
0,0(<0) تباعد 2.8(<>7-,1.8-( 10,12 
0,0(<>31( تباعد 1.8,2.8(<«5)( 20,20 


شكل توضيحي يبين مقارنة بين الطرق من حيث التكرارات 


مقارنة بين عدد التكرارات لايجاد الحل بطرق مختلفة 


(شکل 2 -6) 
ن .ر : طريقة نيوتن رافسون 
ن٠ر.م‏ 1 : طريقة نيوتن رافسون المعدلة 
فر 2 طيقا يرن راقدون امسا وان شا قیھا تددن کم ٭ ایدید 
لايجاد ل 


شکل یوضح تقارب 
النقطة (2,1) إلى 
الحل في طريقة 
نیوتن رافسن 


شکل یوضح تقارب 
النقطة (2,1) إلى الحل 
في طريقة نيوتن رافسن 
المعدلة 


(شکل 2 -8) 


شکل یوضح تقارب 
الشة (01 ال 
الحل في طريقة 
نيوتن رافسن المعدلة 


: باسحل ق چ 
0 1 2 3 4 الجديدة لإيجاد ۷ 
(شڪل 2 -9) 


الأشكال الثلاثة في شكل واحد مع بيان نقطة البدء ونقطة الحل 


خط يوضح تقارب 
الفط را2 ال 
الحل في طريقة 
نيوتن رافسن المعدلة 
الجديدة لإيجاد لإ 


0 05 1 15 2 25 3 35 4 النقطة (2,1) إلى 
ا 


النقطة (2,1) إلى 
الحل في طريقة 
نیوتن رافسن 


4-2 طريقة الموقع الخاطئ 
False position‏ 
sda ERNE EE‏ 
نظرا لأهمية هذه الطرق › رأيت أن أضيف هذه الطريقةء و التي لا تقل أهمية عن 
الطرق السابقة » وللمزيد من الاتساع في الموضوع و لان التعدد يفتح آفاق الخيار 
امام الدارس من اجل التنوع و الاستفادة الشاملة» كنت ازیدذ ان انذا في شرح ) 


interpolation polynomial‏ agrangeا)‏ لا اعرف ترجمة عربية دقيقة لهاء 
على كل حال لنبدا في طريقة ( الموقع الخاطئ( ) false position‏ ( 


المفهوم النظري 
إذا كانت الدالة معرفة في الفترة ( 2× , 1») و كان 2(>0»)؟*(۴)x1‏ 
من الرسم الموضع يمكن استنتاج أن المثلثين : 

0, (x1) , f(x2) 

X3, X2 <2)‏ 
متشابهان فیکون : 
F(x<2)/{f(x2) — f(x<1)}=(x2 — x1)/ (2 —- x1 )‏ 
و تكون قيمة 3× 
X3=x2 —- f(x2)*{(x2- x1)/f(x<2)-f(<2)}‏ 
و يمكن كتابة الخوارزمية للمعادلة التالية 

y1=x**3 +4*x<**2-10 

مع ملاحظة أن الرمز ** يعني في لغة الفورتران للقوة الثانية مثلا أو الثالثة أو 
غيرها آي الاس 

ا 


start .] ® 
read(x1 ,X2) .2 ® 


do .3 ® 
y1=x<1**3 +4*x1**2-10 ® 
y2=x<2**3+4*x<2**2-10 
X3=x2 - f(x2)*((X2- x1)/f(<2)-(<X2)) 
Y3=x3**2+4*x<3**2-0 
if(y3<0.000001)goto (8) .ٍ4 
1f(y1*y3>0.0)goto (7) .,5 © 


xI=×[آ‎ .6 ® 
KISK2 
goto (3) 
KI ل‎ © 
X252 
goto 3 
write x3 ,J3 .& ® 
stop .9 ® 


و البرنامج بلغة فورتران كالأتي : 


Read(*,*)x1 ,x2 
y1=x<1**3 +4*x1**2-10 ® 
y2=x<2**3+4*x<2**2-10 

X3=x2 — f(<2)*((x2- x1)/f(<2)-(X2)) 


Y3=x<3**21+4*x<3**2-10 


if(abs(y3).1t.0.000001)goto 8 
1f(y1*y3.gt.0.0)goto 
ال‎ 

X3SK2 

goto 3 

x1=x3 

RISK 

goto 3 

write(*,*) x3 ,y3 
stop 

End 


F() 


x2 
ا‎ 
1 
+ 
e 


(شكل 2 -11) 
و بهذا نكون قد انهينا طريقة الموقع الخاطئ 


الباب الثالث 


الحلول العددية للمعادلات التفاضلية العادية 


3 -1 الحلول العددية للمعادلات التفاضلية العادية 


سوف نقوم بحل المعادلات من النوع 
dy/dx = f(x,y)‏ 


أي معادلة خطية من الرتبة الأولى . ولحل هذا النوع من المعادلات هناك 
مجموعة من الطرق منها: 

Euler’s Method طريقة اوأر‎ 2= 3 

نفترض أن الحل (×)۴=ر متصل وقابل للتفاضل وأن نقطة البداية هي 
(0,0×) وأن الحل مطلوبا عند ×=« لذلك يجب تجزئة الفترة إلى ١‏ من 
الأجزاء بغرض سا( على أن تكرن صغيرة) 


W = (x -X)/1 حيث‎ 


ثم نطبق القانون 
yı1=¥yÎ +w * f(xi,yi)‏ 
حیث 0,1,..1-1 = 1 
ملاحظة: 
- يكون الخطأ أقل بقيم أقل لس 
- تزداد الحسابات بازدياد عدد الفترات أي تصغیر W۷‏ 
- كل ل تعتمد على التي قبلها مما ينشر الخطأ في حالة وقوعه. 
امتداد طريقة glذر Extended Euler Method‏ 
بأخذ ثلاثة حدود من متسلسلة تايلور التي حصلنا منها على ٤۷‏ نحصل على 
EEM‏ 
yıı=yi +w * f(xi,yi) + w/2! * f `(xi,yi)‏ 
حیث 0,1,..1-1 = 1 
3 -3 طريقة اويلر الكأر ددا More Extended Euler‏ 
Method‏ 
بأخذ الحد الرابع من متسلسلة تايلور نحصJ‏ eلJ” MEEM‏ 
yıı=yi +w * f(xiyyÙ + w/2! * f `(xi,yÙ + w°/3! * f `` (xiyyi)‏ 


Modified Euler Method ةdدعnلا طريقة اویلر‎ 4- 3 


في هذه الطريقة نستغل EM‏ لإیجاد 1+1'ر (٣٥اc ٣٥di‏ ۴) ثم نستخدمها 
اسب ا ا 
avg= ( f(xiyyi) + f(xi+1,y` 1+1) )/2‏ 
نحسب من جدıيد (Corrector) yi+1‏ 
yi+1 = yi + w* avg‏ 


Runge-Kutta Method ligك‎ چiر طريق‎ 5- 3 


نتميز بدقتها وجودتها 
Yi+1= yi + (k1+2k2+2k3+K46‏ 


K1 = w* f(xi,yD 

K2= w* f(xi+w/2,yi+kK1/2) 
K3 = w* f(xi+w/2,yi+K2/2) 
K4 = w* f(xi+w,yi+kK3) 


وبأخذ السؤال الذي يقول أوجد حل مسألة القيم الذاتية : 
dy/dx = xy ; x0=0 , y0 =1‏ 
وذلك عند 1= × علما بأن 0.1=» (قارن بقیم (2/×-)۸ء) 
سنقوم بحله بكل الطرق المذكورة آنفا 
حیث ۷×-=(ل,»)؟ 
اويلر 
Yi+1 = yi - w * xi*yi‏ 
بتقاشل الذالة f‏ وزيادة حذ لايق تحصل غل تداك لويد 
=-xy ¥ y  =-xy +01) = -xy ` -y‏ ¥ 
بالتعويض عن قيمة ل 
Y` = -x(-xy)-y = x y-y= y(x-1)‏ 
الك سل ع مك ور 
Yi+1 = yi - w * xi*yi + w*yi*(xi-1)/2‏ 
بتفاضل ۴ مرة أخرى وزيادة حد نحصل على اويلر الاكثر امتداد 
y= y(x-1) ‡ y ` =y(2x)+(x*-1)y`‏ 
بالتعويض عن قيمة ل 


Y I+1=yi -w*xi*yi 
Avg= ( -xi*yi — xi+1*y`i+1)/2 
Yi+1 = yi + w* avg 


K1= او‎ 


Yir1= + ( COE 


فة نكرر العمل عدد 1 من المرات بعدد الفترات المأخوذة 1 


new‏ : قیمة ر عند بج التالية 


7: قيمة الدالة ( التفاضل) عند ر ,ع 


6. ادخل عدد الفترات ۸ وقيمة لر الابتدائيات 
19. احسب k1=w*F‏ 

F=-)x+w/2(*)y+k1/2( احسب‎ .0 

k2=w*F |حسب‎ .21 

F=-)»+w/2(*)y+٤2/2( احسب‎ .22 

k3=w*F احسب‎ 3 

F=-(x+w)*(y+k3) احسب‎ .24 

k4=w*F |حسب‎ .25 
ynew=y+(k1+2*k2+2*k3+k4)/6 احسب‎ .26 


xıew=x+W |>ب‎ .27 


المخطط الانسيابي 
Flow Chart‏ 


1 1n | 
W—-1/n 


F=-(x+w/2)*(y+K1/2) 
K9I—ut EH 


F=-(x+w/2)*(y+K2/2) 
K3=w*F 


F=-(x+w)*(y+k3) 
K4=w*F 


Ynew=y+(k1+2*k2+2*k3+K406 
Xnew=x+w 


End 
)1- 3 (شكل‎ 
: ملاحظة‎ 
البرنامج الاول هو البرنامج الرئيسي الذي يستخدم طريقة رنج كوتا‎ - 
”=0.1 بالعرض المطلوب في السؤال‎ 

- البرنامج الثاني جعلت فيه W‏ متغيرة بحيث تم أخذ فترات مختلفة 
- البرنامج الثالث استخدمت طريقة اويلر وإمتداداتها 
- البرنامج الرابع طريقة اويلر المعدلة 
- من الممكن طبعا ادماج الطرق جميعا في برنامج واحد 


التعليق 
- طريقة رنج كوتا هي أدق الطرق 
- اويلر المعدلة اعطت نتائج ادق من اويلر الاخريات 
- أما اويلر وامتداداتها فان الاكثر امتدادا كانت أدق ثم امتداد اويلر وأخيرا 
اویلر 


- بزيادة عدد التكرار أي بتصغير ۷ تکون النتائج دق 


0۰6055 => 


(شكل 3 -2) 
العلاقة بين مقدار س وقيمة ر الناتجة حيث يتضح التقارب نحو الحل الصحيح 
بتقليل الفترات 


(شكل 3 -3) 


المحور السيني يمتل قيمة × من 0 إلى 1 بزيادة قيمة W‏ والمحور الصادي هو قيمة 
¥ 
هذا الشكل ييين الطرق لثلاثة. 


3 


منحنى الدالة الأصلي مع القيم الناتجة من طريقة رنج كوتا واويلر المعدلة 


3 -6 طريقة الرمي 


4 


مقدمة 

تستخدم طريقة الرمي لإيجاد حل المعادلات التفاضلية من الرتبة الثانية , ولذلك 
يجب علينا معرفة ثابتين وعادة ما تكون الثوابت إما الدالة ومشتقتها الأولى 
عند نقطة البداية (مسألة القيمة الابتدائية) أو الدالة ومشتقتها عند نقطتين 
مختلفتين (مسألة القيم الحدية) كمسألتنا هذه. 

ان ا ی 

- نفترض الشرط الناقص من شروط مسائل القيم الابتدائية وهو هنا تفاضل 
الدالة عند 1=× 

اوجد الحل بالافتراضات الجديدة (3=×) وقارنه مع الشرط المعطى عند 
- أعد تغيير القيم الابتدائية أي الفرض الذي فرضته حتى نحصل على 
المطلوب وهو (في هذه المسألة 3(=10.0179)ر ) 


کے پک پک پک پک ہے کا ےک ا کے کک ہک پک پیک پیر کے کک ےک کے پک ہک پک ےک ہے کے کک پک ی پک پک پا پیک ہے کے کک ا ا 


وبأخذ السؤال الذي يقول أوجد حل مسألة القيم الذاتية : 
D”y/dx =y ; y(1) =1.1752 , y¥(3)-9‏ 
أولا: نأخذ قيمة افتراضية لر عند 1=× ثم نستخدم طريقة رنج كوتا لإيجاد قيمة ‏ 
نڌ x=3‏ 5 
ثانيا: نفترض قيمة أخرى لر عند 1=× ونعيد استخدام طريقة رنج كوتا لنحصل 
عل ¥ عند 3= , 
الثا نستخدم التوليد الخطي لإيجاد قيمة و أخرى كالتالي: 
VY‏ ۳ 
RI1 g1 g82-g1‏ 
R2 82‏ 
D xo[2]‏ 


Y=y0+pPVy0 
P=(xp-x0)/h 


با o‏ ام 0 د اننا 


H=r2-r1 

VY=92-g1 

P=(d-r1)/(r2-r1) 

۷= ا‎ (d-r1)/(r2-r1)*(92-91) 


ي الخ ف زمية والبر 


tstart‏ : قيمة × ( الابتدائية) 
خاx»sta:‏ قيمة بر الابتدائية 
٥غ‏ متغير يأخذ قيمة ج ويتزايد بمقدار س 


i: R1, r2 
۽ بر ا چ فرعي ل‎ Derives () 


كنا قد تناولنا في المقالات الى 
واحد 2 أن ۰ من المرتبة 


4 -1 طريقة جاكوبي لحل مسائل القيم الذاتية 


مسائل القيم الذاتية تتلخص في ايجاد قيم ۸ وهي القيم الذاتية و × وهي المتجهات 
لذ تتتم الك رة حكر 

Ax= X\Ix 

a1 al2 a13 x1 21 0 0 x1 


a21 a22 a23 |*| x2 |=) 0 22 0 |*| ×2 
a31 a32 a33 X3 0 0 23 X3 


حيث ۸ مصفوفة مربعة متجانسة 
طريقه جاكوبي : 
بهذه الطريقة نريد الوصول إلى المعادل 


Bx= ^X 
حیث ۸ و 8 مصفوفتان فطریتان‎ 
bl1 0 0 
B= 0 b22 0 
0 0 ۵33 
أبعاد < ينة‎ YÈ َد ذات‎ . E 
وطريقة جاكوبي تقوم بدوران المصفوفة بزاوية © في مستوى آخر‎ 
X=Ix 


حيث 1 مصفوفة الدوران ولأنها في ثلاثة ابعاد توجد منها ثلاثة أنواع 


cos -sinO 0 
T1=) sin cos® 0 
0 0 1 
cosټ‎ sinê 
T2= 0 1 0 
sinê 0 cosê 
0 0 
13= 0 cos -sinO 
0 sin cosO 


بما أن 


Ax=\x 
Atx=\ Tx 
۲ بما أن ضرب المصفوفة في معكوسها =1 لذلك نضرب الطرفين في معكوس‎ 
TATx=\T Tx=\x 
ولكي نجعل 747 7 مصفوفة فطرية 8 نختار الزاوية المناسبة أي أن يكون الحد‎ 
الغير الطري في المفوفة الناتجة من الضرب يساوي الصفر‎ 
al2 (cos2 0-sin2 0)+cos Osin O0 (a22-a11 ) 
۴ 
tan (20)=2*a12/ (a11-a22) 
O=tan-1 (2*a1l2/ (a11-a22) ) /2...... (1) 
وللحصول على المتجهات الذاتية نقوم بضرب قات رر ات ال تة من‎ 
ية ى الحضرل عل اق ية‎ 
V=T1*T2*...Tn 


حيث ۸ عدد الدورانات. 


1. نأخذ أكبر عنصر غير قطري في المصفوفة وهو هنا 1 - 
2. مما سبق نأخذ إما 11 أو 13 لأن أكبر عنصر مكرر في مكانين 
3. نحسب فيمة الزاوية من المعادلة (1) وتساوي 
0=tan-1 (-2/ (2-2) =-1/2‏ 
4. نعين المصفوفة الناتجة من الضرب 11۸1 وتصبح هي ۸ 
نضرب مصفوفة الوحدة في ١‏ المرة الأولى ثم الناتج في ١‏ المستخدمة 
أخيرا. وذلك للحصول على المتجهات الذاتية. 
© ت ی ھا کی ی کے ی کی 


:Biger)‏ تحدد أکبر عن 
e‏ تساوي ا 


کے کد الفوز انات 


المخطط الانسيابي 
Flow Chart‏ 


Taoatha—PT/ Taoatha—PT/ 


T44—T( Ua A _- 


ttii=ttjj=tii=t]]=cos(te 
tha) 


A=tt*A*t 


a12&a1 


238.923 


#include<stdio. h> 
#include<conio. h> 
#include<stdlib. h> 


#include<math. h> 
#include<iostream. h> 
#define max 3 


ual (float s[][3], float c[][31]); 
input_mat (float a[max] [max]); 
vei mul1_mat (float aa[] [max], float 
b[] [max], float c[] [max]); 
void write _mat (float a[] [max]); 
void biger (float a[] [max]); 
void tandinvers (float t[][3], float tt[] [3], float 
a[] [max]); 
FILE *in, *out; 
float 
d [max] [max] , a [max] [max] , c [max] [max],tt[3] [3], t€ [3] 


11 "r" ) 


01 open ("jacwrite. fil", "w+"); 

fprintf (out, "The Array of A =\n"); 

a mat (a); 

fprintf (out, "\n");7 

while( (fabs (a[0] [1])>0.0001) I| 

(fabs (a[0][2])>0.0001) || (fabs(a[1] [2])>0.0001) 
{ k=k+1l; 

biger(a ;( 24 find biger element of A 
fprintf (out, "\nCycle no. %d & The biger element 
$5.4f\n",k,a[il[j1); 

tandinvers (t, tt,a ;( // find Angle , T 


fprintf (out, "T = \n"); 
Ww e mat (t); 
mul1_mat (tt,a,d);  // multiplicat invers T in 


mul1l_mat (d,t,a); // multiplicat A in T 
fpPrintf (out, "A = \n"); 

write _ mat (a); 

mull _mat (v,t,c); //clculat Eigenvectors 


} 
۶( چ ane‏ اوا و 


E 45. 4£ ",v[iIl[i1); 
J 

fclose (out); 

return (0) ۶ 


// function to find angle , T and invers T 


void dinvers (float t[][3], float tt[] [3], float 
a[ ] [max]) 
{ float x,z,y,tunit [3] [3]={1,0,0,0,1,0,0,0,1}; 
1if(a[i] [i]==a[j][j1]) if (a[i] [i1>=0) x=3. 1416/4; 
else x=-3.1416/4; 
ee x= (atan (2*a[i][j]/(a[i] [i]-a[j][j1])))/2; 
y=cos (x); Eee SF 
equal (t, tunit); // t= {1,0,0,0,1,0,0,0,1} 
equal (tt, tunit); 
// select T1 , T2 or T3 consider biger element 
4 
ج‎ =0;i<n; i++) 
or (j=0; j<n; j++) 
fscanf (in, "$f", éa[il[j1l); 
} 
void biger (float a[] [max]) 
{ 


if (fabs (a[0] [1] ) >=fabs (a[01] [2 (([ 
if (fabs (a[0] [1])>=fabs(a[1][2])) {i=0;j=1;} 
else {1=1; =2; [ 
else if (fabs (a[0] [2]) >=fabs (a[1] [21])) 
11=0; =2; [ 
else {i=1; =2; 


J} 
void mull _ mat (float aa [] [max], float 
b[] [max], float c[] [max] ) 
{int k; 
for (i=0; e 


oy r 
c[i][j]=c[i] [j] +aa[i] [k]*b[k] [j] 2 


7 


void write_mat (float c[] [max]) 


{ 


for (i=0;i<n; i++) 


{ 


for(j=0; j<n; j++) 
printf (out, "%5.4f\t",c[i][j1); 
fprintf (out, "\n"); 


0.0000 
1. 0000- 


The Array of A= 


1. 0000- 


1 & The biger element : -1.0000 


0.7071 
2. 0000 


0. 0000 


0. 0000 
1. 0000 


2 & The biger element : 0.7071 


0.4597= 
0. 0000 
0.8881 


0. 0000- 


0.0000 
1. 0000 
0 . 0 0 00 


0.3251- 
1 . 0000 


3 & The biger element : -0.6280 


0.0000 
0.5241- 
0.8517 


0.0000 
0.8517 
0.5241 


2. 0000 
1. 0000- 


0.4597 
4 5 

3. 0 
00.3251- 
0. 0000- 


J} 


0.1695 


0.0020 


. 4 
. 0000 


. 0204 


0 
0 
٤ 


0.2768- 
0.6136 
0. 0000 


3.30 


0.2768- 


0.1704 


Cycle no. 4 & The biger element 


0.0991 
0.9951 
0 : 0000 


0. 0000- 
0.5860 
0.0169 


5 & The biger element 


1207 


. 0000 


. 97 


.0168 


. 998 


O O O 


HH O O 


0. 0000 
1. 0000 
0. 0000 


0.5860 
0.0168 


6 & The biger element 


cle no. 7 & The biger element 


Oo OO 


0.0004 
۳ 


0.0000 
0.9999 
0,0119= 


2 
585 


0.0007- 
1 . 0 0 O00 
0.0000 
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.0991- 


. 0000 


coum oo 
د‎ © u اا‎ 


T = 


J0=3. 42 
V= 0.5000 -0.7071 0.5000 


J1=0.5858 
V= 0.5000 0.7071 0.500 
J2=2. 0000 
V= -0.7071 -0.0000 0.70771 
الاستنتاج والتعليق‎ 
تم استخدام عدة برامج فرعية لقراءة وطباعة المصفوفة ولحساب الضرب‎ = 
وغیره‎ 


- يختلف عدد الدورانات باختلاف الدقة المختارة للخطاً المسموح به فكلما 
زادت الدقة كلما زاد عدد الدوران فالعلاقة طردية (عند 0.0001 عدد 7 
دورانات) 

- نلاحظ في حالة تساوي قيمتين كأكبر قيمة فإن الاختلاف في أخذ 1 يؤدي 
إلى نفس الحل ولكن ليس بنفس الترتيب . 

- استخدمت شرطا خاصا لإيجاد قيمة الزاوية عندما يكون المقام بصفر 

- لاحظ أن من موقع آكبر عنصر نستطیع تحدید ۲ فلو كان موقع أكبر 
عنصر =ز,ا=آفإننا نضع all=cos, akk=cos , alk=-sin , akl=s1¬‏ 

- استخدمت مصفوفة الوحدة لوضعها في مصفوفة الدوران قبل تحميلها كما 


. 


سبق . 


4 -2 طر يقة 1ظ جر انج للتو ليد Lagrange’s Interpolating Polynomial‏ 


ر ا ر غ جه ت ی جو قرات 
الحدود وكذلك لايجاد صيغة كثيرة الحدود . وهي تصلح لأي ند ھن 
النقاط. وتعتمد درجة كثيرة الحدود الناتجة على عدد النقاط المستخدمة. 
درجة كثيرة الحدود = ۸-1" حيث 1 عدد النقاط 
فلو كانت عدد النقاط 3 فستكون كثيرة الحدود من الدرجة الثانية وهكذا 
وتكون عدد معاملات لاجرانج ] تعتمد على عدد النقاط وتكون قيمة ا 
على الصورة: 

Lis _ = DC = x2... = Xi =D = i + D.x... (= r) I 


(Xi = Xi = X2)......(%i = Xi = 1)(Xi = Xi + 1.......(Xi = Xn) 
ححاصل ضرب حاصل طرح قيمة ×من كل النقاط عدا النقطة الحالية‎ 
حاصل ضرب حاصل طرح النقطة الحالية من كل النقاط عدا نفسها‎ 
: تم لإيجاد المعادلة النهائية‎ 
Y=L(x) = JLi * yi 
وهن المااكظ انه لا يمك بصهولة أى قد بر لها يريا خاصة عنما‎ 
تكثر عدد النقاط ويزداد صعوبة عندما تكون قيم النقاط بالأرقام العشرية.‎ 
© لذلك تحل باستخدام الحاسوب. وسنبرمج طريقة لاجرانج باستخدام لغة‎ 


قائمة بالمتغيرات المستخدمة في البرنامج 

+ کرات رار 

د التقاط المنخلة 

صر : النقطة المراد إيجاد قيمة ر عندها 

فر ,4× : النقاط المدخلة 

)1( حاصل ضرب طرح قيمة ص من النقاط أي البسط في القانون‎ :٤ 


ل زر الثالى وي . 
هل انتهت النقاط )ص 


المخطط الانسيابي 
Flow Chart‏ 


No 


D=D*(xp-xi) 
CC *(wi-wi) 


# include<stdio. 


# include<conio. h> 
void main () 
{int i; jy} 

float L[9], sum 

NSO 
printf ("Enter n value\n "); 
scanf ("%d" , en); 
printf("\n Enter x and y \n "); 
for (i=0;i<n; i++) 

scanf ("%f Sf",&x[i],&y[i]); 

printf ("\n Enter point of x "); 
scanf ( 2 ", &xXP) 
print 2 hê method\n\n") 7; 

printf ("\n 5 

for ( 1<n; 31++ 1 printf ("%4.2f EEE 

for =0; زو ا‎ PFiRntEF("%$4.2F\E", III; 
printf (" xp= %4.2f\n", xp); 
E ر‎ 
n"); 

for (i=0;i<n; i++) 

for (j=0; j<n; j++) 
C=C* (xp~x[j]) ; 
D=D* (x[i]-—x[j]); 
EN e 
printf("\n L(%d)= %f ",i,L[i]); 
j} 

£("\n\ 

n"); 

he EÊ ("\nl=L1*y1+12*y2+....\n"); 

printf A 7 
RT} 

printf ("L= $f ",sum); 


ء 


1 | * 
1 N 
1 ملار‎ 


Lagrange method 


x 0.00 1.00 2.00 3.00 4.00 
y 0.00 2.00 8.00 18.00 32.00 xp= 


1(3)= -10. 000000 


I= 50.000000 


و sin(x)‏ تکون النتائج 


Lagrange method 


E E 90.00 


0.00 0.50 0.8660254038 1.00 xp= 


1L(1)= 0.740741 
1(2)= 0 


L=L1*y1+12*y2+.... 


1= 0.641738 


1 (0)= -0.049383 
1(1)= 0.740741 
1L (2)= -0.061728 
1(3)(= 0 


Lagrange method 


40.00 


1(0(= -.1 
1L(1)= 0.888889 


1(2) =  . 2222 


E 


a e 


30.00 60.00 90.00 
0.50 0.8660254038 1.00 xpD= 40.00 


* 


E 0.555556 
1)2(= - - 1 


L=L1*y1+12*y2+.... 


كه . ففي 2- ع كانت ص× قريبة من النقطة الث 
كانت ط× قريية من النقطة الأولى في 2 -ل زادت ا[ فهي 


Dimension x(50),f950) 
Write (*,*)'enter value of polynomial’ 
Read (*,*) n 
Write (*,*)'enter point of x that you find solution’ 
Read(*,*)xO 
Do 10,i=1,n+1 
Write (*,*)'x (iD), fi, )=' 
Read(*,*)xD,fO) 
10 Continue 
Do 20,i=1,n+1 
t=1 
Do 30,j=1.n+1 
If(i.ne.j)then 
_ T=(t*(x0-x(j)))/(x(i)-x(j)) 
Else 
Endif 
30 Continue 
Px=px+t*(f@ 
20 Continue 
Write(*,*)xp 
Stop 


3-4 حذف جاوس: طريقة التعويض الخلفي 


قبل الدخول في التطبيق البرمجي › أرى من الأهمية بمكان طرح المفهوم النظري 
لطريقة جاوس أو طريقة ( التعويض الخلفي) . 


الأساسيات النظرية 


تقوم طريقة جاوس على تحويل المعدلات إلى نظام المصفوفات و التعامل معها 
على هذا الأساس كما الآتي : 


A1x1 + a2x2 + a3x3 +....... anxn = c1 


B1x1 + b2x2 + b3x3 + ......... bnxn= c2 


M8 a. I bê ba i a را‎ 


X1 


“Q 


مصفوفة المرافقا 
بمقدار واحد , لتصڊ 


al,1 a1l,2 .... aln c1 


F2,1=b2,1/a1,1 
b2,1= b2,1 - f2,1*a1,1 
21*a1,2 


b2,2 =b2,2 - f 
b2,3=b3,2 - f1,2*a1,3 


الذي قمنا به هو الآتي : 


4)50,50( عرف المصفوفات (50) × (۴)50 (50 ,50) ط‎ 
read n 
do 
from I = 1 to n read elements of matrix a 
do 
m=n*(n-1) 
do 
from I=2 to m 
from k=1 to i-1 
if (a (k,k) #0 yes f(I,k)=a(1,k) / a(k,k) 
from j = 1 to n +1 
a(1,j)=a(1,j) - f(1,k)*a(k,j) 


no 


from |1 =I — 1 to I —1 
from j=1 to n+1 
b( ,j)= a1 ,j) 
a(l,j)=a(1+1,j) 
al1l+1,j)=b (1,j) 
do .9 
from I = n to 1 step —1 
from j =I+1 ton 
s1= a(l, n +1) 
x(i) = (s1 —- s2)/ a( 1,1) 
do .10 


from I = 1ton 


write x(i) 


(n +1 


البرنامج بلغة(الفورتران ( FORTRAN language‏ 
(50)× , (۴)50 , ( 50 ,50) طاو 


Dimension a(50,50) 


Write(*,*)' enter numbers of your elements (n)' 
Read (*,*)n 
Do 10,i=1,n 
Do 20 , j=1 n+1 

Write(*,*)' a ONY =" 
Read(*,*)a(1,j) 

20 continue 

10 continue 
m=n*(n-1)/2 
do 30,i=2,m 
do 40,k=1,n-1 
if(a(k,k).ne.0)then 
do 40,j=1.n+1 


a(i.j)=a(1,j)-fI,k)*a(k,j) 
40 continue 
else 


do 50,L=i-1,i-1 


do 60,j=1,n+1 

b(1,j)=a(1,j) 

a)1,ز(=21+1,ز(‎ 

a)1+1,j(=b)1,ز(‎ 

60 continue 
50 continue 

do 70,i=n,1,-1 

do 80,j=i+1,n 

(j) 


80 continue 


S2=s2 + a(1,j)* 


s1=a(l,n+1) 

x(i)=(s1 — s2)/a(l,i) 

70 continue 

do 90,i=1,n 

write(*,*)' x ('jl,' ) = 'yx(D 
90 continue 


stop 


end 


A(Lj) =a, j )-f (1 k)*a(k, j ) 


a(1,j)=a(1+ 1 , j) 


al + 1, j(=b)1, ز‎ ( 


أ ”1° *1 Ah‏ 
ډلیء 
چ 


نكون متغير يحمل الاسم 2ء الذي يقوم بقسمة | 
العنصر الذي رتبته ز و كأننا قمنا بالاتي'" 


۵٤‏ ا 1 ! ف S2 1 ۳ a 8 i ٤ ۴ ١‏ و جه ١‏ 2 الاستمرا 
0 چ ٠ e‏ 
1 سو _ 


سل | اسه 
+ 
٤ a . TT,‏ 
TT 1‏ 2 ا TIT‏ 
EE‏ أللےعه لص 


“FITTING METHOD 


FITTING METHOD cii andl 1-5 


عندما تعطى بيانات غير دقيقة نوعا ما أي تقريبية كالقياسات المعملية 
والهتدسية فان هذه البيانات ثر تبط بأخطاء مضدذرها القياس أو الاسان أو : 
لذلك نلائم المنحنى الناتج من هذه البيانات حيث أن هذه البيانات تتبع معادلة معينة 
القن دنام نفک کک من ق رقن انی ارآ رای ن 
سضر البيانات نها هل تمل المسقلة ذالة تر ية أو اسية أو ١‏ وهكذا وبعذ 
تيد الذالة يجب أن نجد اجرد ملمة ۴۲ )ئ86 وهي تعني أن نجعل الأخطاء 
أقل مايمكن فنأخذ مجمو ع مربعات الأخطاء 


S= 3 
i1=0 


8 = Pp(x)-yi 
Pp (%) = ag+aıxX +....... +aqX" =), aj£8j(X) 


ا ي غر ا م وات كف ارات ار 
as‏ 
وتكون $ آقل ما يمکن عندما ٠‏ چ 


وسندخل الآن في متال عملي على ذلك 

الجدول أسفله يعطي بيانات عن مدى استهلاك الماء بإحدى البلدان وذلك 

ببلايين الجالونات في اليوم 

- استعمل الانكفاء الأسى لاستهلاك الماء بدلالة الزمن 

- استعمل ما حصلت عليه لحساب استهلاك الماء في سنة 1975 وقارن بما 
کان متوقعا وهو 449.7 


1970 1960 1950 1940 1930 السنة 
411.2 32.9 27 || 136.43 5 | الاستهلاك 
(جدول 5 -1) 


وسنطبق ما سبق على المعادلة التالية (الاسية) 
Y=ab*‏ 
ولكي تصبح خطية نأخذ لوغارثم الطرفين 
Ln (y) = In ( ab®)‏ 
Ln(y) = In(a) + x In(b)‏ 
Z=۸A+Bx‏ 
حيیث B=1n(b)‏ و A=1n(a)‏ ڍ Z =1n(y)‏ 


$= 3 _ J (A+B»-2 
i=0 i0 


dS n 
— =2 (A+Bx- ك‎ 
3 +Bx-2)()=0 


dS n 
— =2 )A+Bx- =0 
3 +Bx-2)(9) 
المعادلات العمودية‎ 
nA ۴ B}, Xi = 3 Zi 


A), Xt B}), x = 3 ZiXi 


7Z =[۸)y( وبارجاع‎ 
nA B}), Xi = 3 In(y; 
AJ xi+ BY x7 = J x;ln(yi) 


n 2 
۸ = غ ت‎ 7 (1) 
JJ xi E 
JJ In(yi) E 
> /۸ = [}, xi, In(yD- }D, xi }, xiln(yDY 
JJ xiln(yî) 3 
Al) 
n JJ In(yi) 
B= /۸ = [n}, xiln(yD- DJ In(yD D, xil/ ^ ....... 


3 JJ xiln(y) 


a=e* , b=e% 
بالعودة للمعادلة الرئيسية والتعويض بقيم 4 و طط‎ 


TP ANGERS (4) 

وطبقا للبيانات المعطاة والتي تبين مدى استهلاك الماء بإحدى البلدان ببلايين 
الجالونات. 

1970 1960 1950 1940 1930 السنة 


411.2 32.9 2027 1| 136.43 5 1 الاستهلاك 


(جدول 5 -2) 
ولو أردنا إيجاد الحل يدويا فنكمل الجدول ونختصر السنوات كالتالي: 


x |y Ln(y) | X X*In(y) X*In(y) 
30 110.5 4.7 900 141.15 4234.51 
40 136.43 | 4.9 1600 196.63 7865.30 
50 202.7 5.3 2500 265.59 13279.32 
60 32.9 5.7 3600 346.64 20798.43 


70 411.2 6.0 4900 421.4 2943.49 
250 1183.73 26.6 | 13500 /ً|nO1C1"1 1.5 75671.05 
)3- 5 (جدول‎ 


بالتعويض عن قيم المجاميع 
5000 = 250(7) -5*13500 = ۸ 

A= (26.6*13500-250*1371.35(/5000 =5 

B = (5*1371.35-250*26.6)(/5000 =0. 041 

a=مe‎ 25 = 25,79 


Y=ab* 

Y=25.79*(1.04)* 

J (75)(=25.79*(1.04)(” =501.8 

هذه الأرقام العشرية ليست دقيقة بل مقربة لرقم أو رقمين ولكن النتيجة النهائية 
E‏ 

لرل ك اللاتة 
70 60 50 40 30 
4.46 297.31 209.72 147.94 104.36 


(جدول 5 -4) 
قائمة بالمتغيرات المستخدمة في الخوارزمية والبرنامج 


1 متغير للتكرار 

N‏ عدد النقاط المدخلة 

ذس × : قيم د(السنة) ومر(الاستهلاك) المعطاة 

عد السخة المر اد ايجك قيمة الأستهلاف عندها 

: قيمة الاستهلاك الناتج عند ب 

2غ 1عل: قيمة المحدد العام الناتج من المعادلات السابقة...... المعادلة (1) 


×صای: مجمو ع قیم × 


64. احسب مجga‏ ع (sumxy=sumxy+xi*1n(yi)) xi*1(yi)‏ 
6. هل >1 إذا كان نعم ارجع إلى 4 
7. احسب قيمة 4٤امل‏ من (1) 4 من المعادلة ( 
. ادخل قيمة × المراد ايجاد ل عندها 
لفت ق ون ا ۲ 
0. النهاية 


المخطط الانسيابي 
Flow Chart‏ 


SUMX=sSUMX+x 
sumy=sumy+y 


SUMXX=SUNnIXX+X* 


sumxy=sumxy+x*In(y 


T—-Tı1 


4 
: © 
س 


#include<stdio. h> 
#include<math. h> 

void main () 

{ 

int i,n, sumx 
double a,b, sumy, y, sumxy, 


x, sumxx, xi1[5];‏ ر 
elta,yi[5];‏ 


printf ("Enter Points Number: "); 
printf ("Enter Values of x and y "); 
sumx=0; sumy=0; sumxy=0; sumxx=0; 


for (i=0;i<n; i++) 

{ 

fscanf ("%d", &xi[i]); 
scanf ("%lf",&yi[i1]); 


(yi i1);‏ و ر 


e 1 


printf ("Enter value of x"); 

scanf ("%d",&x); 

pow (b,x);‏ ر 

// writing 

printf ("a= $ , b= $1f\n\n",a,b); 
printf(" y = a * b ^ x\n\n"); 
printf ("x : Ji<n;i++) printf("%d\t",xi[i] ); 
prin (ay "J; for(i=0;i<n;i++)printf ("%5.21lf 
",a*pow(b,xi[i1])); 
printf("\n\n y = %5.2lf * %5.21f^%d = %1f\n\n",a,b,x,y); 
printf ("When x=%d then y=%5.3lf",x,y); 

J 


® 

1 

| oo 

m0 H 
N 

1 


البرنامج: في صورته وهو يقرأ من ملف ويطبع إلى ملف 
#include<stdio. h>‏ 
#include<math. h>‏ 
void main ()‏ 
{ 
int i,n, sumx, x, sumxx, xi[5];‏ 
double a,b, sumy, y, sumxy, delta, yi[5];‏ 
FILE *stream;‏ 
input values‏ // 
stream = fopen ("fitread.FIL", "r");‏ 
fscanf (stream, "$d", en);‏ 
for(i=0;i<n; i++) scanf‏ 
for(i=0;i<n; i++)‏ 
fclose (stream) ;‏ 


sumy=sumy+yi [i]; 


sumxx+xi [i1] *xi [i]; 
y=sumxy+xi [i] *yi [i]; 


delta = RSE SUE 


prin ee of × me 
scanf ("$d", &x) 7 


pow (b,x); 

rIting 

= fopen ("fitwrite.FIL", "w+"); 

tf (stream, "a= %$lf , b= $%$l1f\n\n",a,b); 

tf (stream, " y = a * b ^ x\nl\n"); 

و ر" am,‏ 

Foz (=0, i<n; i++) fprintf (stream, "$d\t",xi[i] ); 
fpri (stream, "\ny :"); 

for(i=0;i<n;i++) fprintf (stream, "%55.2lf ",a*pow(b,xi[i] )),; 
fprintf (stream, "\n\n J = %5.2lf * %$5.21f^%d = 
$1f\n\n",a,b,x,y); 

fprintf (stream, "When x=%d then y=%5.3lf",x,y); 
fclose (stream) ; 


j 


الإادخالاتثت 
Enter Points Number: 5‏ 
En er Value of x 30 40 50 60 70‏ 
r Value of y 110.5 136.43 202.7 322.9 411.2‏ :1 
Enter value of x 75‏ 


a= 36.633592 , b= 1.03 


70 60 50 40 30 :2 ی 
yJ :104.36 147.94 209.72 297.31 421.46‏ 


yJ = 36.63 * 1.04^75 = 501.87 


When x=75 then y=501.804 


النة 
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استهلاف 


لاء 


اأسة ةة 


النقاط المعطاة مع منحنى الملائمة 
(شکل5 -3) 


في الختام لا يسعنا الا أن نسأل 
الله عز وجل تحقيق الفائدة من هذا 
الكتاب ونتمنى آنه لقي استحسانكم . 
ویسرنا أن نتلقی آراءکم حول فحوی 


الكتاب على عنوان البريد الالكتروني : 


tkne@tkne.net 
ولكم جزيل الشكر‎ 


وادارة موقع النقنية 


